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Abstract 

By making use of the finite-temperature real-time static potential that was introduced and 
computed to leading non-trivial order in Hard Thermal Loop resummed perturbation the- 
ory in recent work, and solving numerically a Schrodinger-type equation, we estimate the 
quarkonium (in practice, bottomonium) contribution to the spectral function of the electro- 
magnetic current in hot QCD. The spectral function shows a single resonance peak which 
becomes wider and then disappears as the temperature is increased beyond 450 MeV or so. 
This behaviour can be compared with recently attempted lattice reconstructions of the same 
quantity, based on the "maximum entropy method", which generically show several peaks. 
We also specify the dependence of our results on the spatial momentum of the electromagnetic 
current, as well as on the baryon chemical potential characterising the hot QCD plasma. 
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1. Introduction 



It was suggested long ago that the properties of heavy quarkonium may be very sensitive to 
the deconfinement transition that takes place in thermal QCD, in spite of the fact that the 
deconfinement temperature is much below the heavy quark mass PQ. Consequently, heavy 
quarkonium has become one of the classic probes for quark-gluon plasma formation in heavy 
ion collision experiments (for an extensive review, see ref. [2]). 

In order to understand the physics involved, let us start by recalling that the way in which 
the properties of thermally produced heavy quarkonium can be observed, is through its decay 
into a virtual photon, which then produces a lepton-antilepton pair [3]. Leptons do not feel 
strong interactions, and escape the thermal system. Measuring their energy spectrum at 
around E ~ 2M, where M is the heavy quark mass, thus yields first-hand information on 
the "in-medium" properties of heavy quarkonium. 

To appreciate why the in-medium properties of heavy quarkonium can change already just 
above the deconfinement transition, it is conventional to consider a non-relativistic potential 
model for determining the thermally modified energy levels of the decaying bound state [4]-[9]. 
Above the deconfinement transition, the colour-electric field responsible for binding the heavy 
quark and antiquark together gets Debye-screened. Once the screening is strong enough, the 
corresponding Schrodinger equation does not possess bound-state solutions any more. It is 
said that quarkonium "melts" at this point, and the assumption is that the quarkonium 
resonance peak should consequently disappear from the dilepton production rate. 

Strictly speaking, though, just estimating the energy levels from a potential model does 
not allow to reconstruct the spectral function (which in turn determines the production rate) . 
In fact, stationary levels would correspond to infinitely narrow peaks in the spectral function, 
irrespective of the value of the binding energy, while the intuitive picture is that a resonance 
peak should dissolve through becoming gradually wider. To conform with this expectation, 
a non-zero width could of course be inserted by hand, as an additional model ingredient. 
However, this would take us further away from a first principles QCD prediction. 

It appears that once the computation is formulated within thermal field theory, there is no 
need to insert anything by hand. Indeed, it has been pointed out recently that by defining 
a static potential through a Schrodinger equation satisfied by a certain heavy quarkonium 
Green's function, and computing it systematically in the weak-coupling expansion (which 
necessitates Hard Thermal Loop resummation), the static potential obtains both a standard 
Debye-screened real part, as well as an imaginary part, originating from the Landau-damping 
of almost static colour fields [10] . The imaginary part of the static potential then leads to a 
finite width for the quarkonium resonance peak in the spectral function. 

In ref. [10] . the consequences deriving from the existence of an imaginary part were ad- 
dressed only semi-quantitatively. It is the purpose of the present note to solve explicitly for 
the spectral function that the static potential computed in ref. [10J leads to. We also compare 
qualitatively with attempted lattice reconstructions of the same quantity. 
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The note is organised as follows. We review the form of the spectral function in the non- 
interacting limit in Sec. [2j Some properties of the static potential derived in ref. [10] are 
discussed in Sec. El The relevant (time-dependent) Schrodinger equation is set up in Sec. HJ 
and solved numerically in Sec. We conclude and compare with literature in Sec. El 



2. Spectral function in the non-interacting limit 

We will consider two related correlators in this paper: 

C>(q°) = l°°Jt J d 3 xe^(^(x) ^(O)) , (2.1) 

where j^{x) =ip {x)^t{){x) is the contribution from a single heavy flavour to the electromag- 
netic current in the Heisenberg picture (the electromagnetic coupling constant and charge 

have been omitted for simplicity, and the metric is assumed to be (H )); as well as the 

spectral function 

p(Q )^\(l-e-^)c>( q °), (2.2) 

where (3 = 1/T, and T is the temperature. The dilepton production rate is directly pro- 
portional to the spectral function [3|. The expectation value in Eq. (|2.ip refers to (...) = 
2 _1 Tr [exp(— H /T)(...)], where Z is the partition function, and H is the QCD Hamiltonian 
operator. We have assumed a notation where the dependence on the spatial momentum q 
is suppressed. A correlator without tilde refers to the situation before taking the Fourier 
transform with respect to time: 

C>(t) = |d 3 xe-^ x (^(t,x)^(0,0)) . (2.3) 

We start by discussing the form of p(q°) in the free theory. Denoting the heavy quark mass 
by M, we concentrate on frequencies around the two-particle threshold, 



w = q° y4M 2 + q 2 , (2.4) 
and will also assume the spatial momentum q to be small, q = |q| <C M. 

2.1. Non-relativistic low-temperature regime in full QCD 

The free quarkonium contribution to the spectral function of the electromagnetic current can 
be extracted, for instance, from refs. [11 } 112 } 19]. Modifications brought in by various lattice 
discretizations have also been addressed [111 \12\ I13j . Here we generalise the continuum 
expression slightly by including a non-zero quark chemical potential, [i. Restricting first to 
the case q = 0, the result is very simple: 



j>M N c ,, 2 /i/~ on/i 4 



P {u) ^= -jZMPeifi - 2) (1 - ^ J - {& + 2)[1- n F (| + M ) - n F (| - M Jj , (2.5) 
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where N c = 3, n-p is the Fermi distribution function, and we have denoted 



M 



(2.6) 



Let us now concentrate on the case of low temperatures, T/(M ± //) <C 1 (parametrically, 
we are interested in temperatures T ~ g 2 M |10j). Then the Fermi distribution functions in 
Eq. (|2.5p are exponentially small. We thus find immediately that the spectral function is 
independent of fi in this limit. 

Restricting furthermore to the non-relativistic regime, \u — 2| <C 1, and considering the 
external momentum q to be small, q <C M, it is easy to include dependence on q. We obtain 



o^2M 3N C M 2 



PV ' 2n V 4M 2 A 4M 2 / 



i r 

2 



1 + 



9 



AM 2 ' M 2 / 



• (2-7) 



2.2. Representation through a Schrodinger equation 



We next demonstrate that the result of Eq. (|2.7p can be reproduced by a certain Schrodinger 
equation. The Schrodinger equation requires the introduction of an intermediate point- 
splitting vector r which will be set to zero at the end of the computation. The relevant 
equation reads (cf. Eq. (2.4) of ref. [10] ) 



id t 



2M 



2\ i 



C>(i,r) = 
6iY c 5 (3) (r) . 



(2.8 



with the initial condition 

C7>(0,r) = -6N c 5^>{r) . (21 

In Eq. (|2.8p we have set for simplicity q = 0, but the center-of-mass kinetic energy q 2 /AM 
can be trivially added to the rest mass 2M. After having solved the equation, the function 
in Eq. (12. 3p is obtained through 

C>(i) = £>(t,0). (2.10) 



(2.11) 



We search for a solution of Eq. (12, 8p with the ansatz 

d 4 P 



C>(i,r) 



(2*)< 



e - W t+ipT f(p(Jip) . 



Eq. (g2D dictates that 



p ' 



leading to the modified ansatz 

C>(t,r) 



I (W e "" pt+ "'^( p ) • 



(27T) 



(2.12) 



(2.13) 
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The initial condition in Eq. (|2.9p can be satisfied provided that J-(p) = — 6iV c . The point- 
splitting can now be trivially removed, cf. Eq. (|2.10p . and a Fourier-transform finally yields 



/oo 
die^C^O) 
-oo 



3iV r M 2 



(2vr) 3 V M 



6(u> - 2)(w - 2)2 , (2.14) 

7T 

where we have used the dimensionless variable in Eq. (|2.6p . The spectral function is given 
by Eq. (j2.2|) : since we are in the non-relativistic limit — 2| <C 1 and at low temperatures 
T <C M, the factor exp(— /3u) ~ exp(— 2M/T) can be neglected, whereby p(w) = C > (u>)/2. 
Replacing furthermore 2M — > 2M + q 2 /AM, yields then directly Eq. (j2.7[) . as promised. 



3. Real-time static potential 

In order to account for interactions, a static potential can be inserted into the Schrodinger 

(2) 

equation. The appropriate object, denoted by V> (t, r), was defined and computed to leading 
non-trivial order in Hard Thermal Loop resummed perturbation in ref. [10] (cf. Eq. (3.17)). 
Reorganizing the result in a way where the symmetry of the integrand under p° «- 
explicit, we rewrite it as 



p° is 



Vf } (i,r) 
5V± 2) (t,r) 



9 2 C F 
4tt 



g 2 c F 



m D + 



exp(— m D r) 



d 3 p 2 



(3.1) 



(2vr) 3 
dp o 



-P 



e -*|p°l* + n B (|p°|)(e- < l^l* - e*^l* 
^2 ~ J^2-)pe(p°iP) + - p)^r(P°,P)] } 



(3.2) 



Here (7f = (iV^ — 1)/2N C , m D is the Debye mass parameter, and we have chosen r = (0, 0, r). 
The r-independent term in Eq. (|3.ip amounts to twice a thermal mass correction for the 
heavy quark. The functions pe,Pt are specified in Appendix A. The Schrodinger equation 
to be solved reads 

idt - ( 2 M-^ + vi 2 \t,r) 



C>(t,r) = 



(3.3) 



with the initial condition in Eq. (|2.9p . and the replacement 2M — > 2M + q 2 /AM for <7 7^ 0. 



3.1. Dynamical scales 

Let us review the time and distance scales that play a role in the solution of Eq. (|3.3p . The 
derivatives in the free part must be of similar magnitudes (after trivially shifting away the 
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constant 2M), implying that 

1 fl\ 2 1 

(3.4) 



t \r J M 

At the same time, they must also be of similar magnitude as the potential. Given that the 
potential is screened, this means 

r J M r r 
Therefore, we obtain 

7 S ?\ ■ (3-6) 

i.e. the time scales relevant for the solution around the resonance peak are much larger 
than the spatial distance scales. Consequently, in order to obtain a formally consistent 
approximation to a fixed order in g, we need to take the limit t S> r in the static potential. 

Even though it has thus become clear that only the limit t S> r of the potential is needed 
at the first non-trivial order in g 2 , we nevertheless discuss in the remainder of this section 
how the infinite-time limit is approached, perhaps learning on the way something about the 
convergence of the weak-coupling expansion. 

3.2. Zero-temperature part 

Let us first compute 8V>\t,r) in the zero-temperature limit. In this case nB(|p°|) — ► and 

PE(P°,P) = PT(P°,P) =7Tsign(p )5((p ) 2 - p 2 ) . (3.7) 

Given that the prefactor in front of pe vanishes on-shell, pe does not contribute in this limit, 
and we simply obtain 

where p = |p|. Even though it is obvious that this contribution vanishes for t — > oo, its precise 
evaluation requires the introduction of an intermediate regulator, because the absolute value 
of the p-integrand grows linearly with p. We can either set t t — ie, with e — > + at the end 
of the computation, or regulate the spatial integration by going to d = 3 — 2e dimensions. In 
the first case the integral can be rewritten as 



Svi 2 \t,r) = ^ j\z (1 - l) J^dpe-^e-^ - e^ rz ^] ; (3.9) 

in the latter case the "convergence factor" e~ pe is replaced by p~ 2e . Either way, the p- 
integral can be carried out (in the former case, dpe~ pe e~ ipx = l/(ix + e); in the latter 
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Figure 1: The part of <5V> (t, r) that remains finite for t — > oo (cf. Sec. 13. 3p . The circles at right 
denote the asymptotic values in this limit. The oscillations visible at large rm D have the frequency 
djpi = m D /-\/3; the corresponding oscillation period in terms of the variable im D is 2n^/3 « 10.9. 



case, / °° dpp 2e e tpx = T(l — 2e)/(ix) 1 2e ), and subsequently, also the z-integral (as long as 
we stay within the light cone). We obtain, for t > r, 

2 + - 1 - — In ■ 



5vi 2 \t,r) = g 2 C F 



4vr 2 i 



XT' 

9 2 ° F ^3> for t » r . (3.10) 



t\ r 2 J t — r 

The result is, thus, purely imaginary, and vanishes rapidly with time. For t^ 1 ~ g 2 r~ l , it 
corresponds parametrically to an effect of order 0(g 8 /r), and should be neglected. 

3.3. Finite-temperature part 

(2) 

Considering then £V> (t, r) at finite temperatures, there are two different types of new struc- 
tures emerging. First of all, there is the term without ne(|p |) in Eq. (|3.2p . This amounts to 
a generalization of the potential in Sec. 13.21 through the introduction of one new dimensionful 
parameter, m D , appearing in the spectral functions. Second, there is the term with 71b(|p°|)- 
This introduces a further new dimensionful parameter, T, and complicates the functional 
dependence further. 

The evaluation of the term without nedp ]) again requires the introduction of a regulator, 
as in Sec. 13.21 The resulting potential has both a real and an imaginary part. However, it 
still decays fast for t^>r; the only difference with respect to Sec. 13.21 is that the decay is not 



6 



purely powerlike any more, but the existence of a new scale leads to oscillations as well. In 
particular, at large r the behaviour is dominated by small p, and then the oscillations take 
place with the familiar plasmon frequency, oj p \ = m D /v3 (cf. Eqs. (jA.llj) . (IA.12j) ). 

On the other hand, the term with nB(|p°|) leads to more dramatic new effects. As is obvious 
from Eq. (|3.2j) . the contribution from this term to the static potential is purely imaginary. 
Also, this part can be evaluated without regularization, since nedp )) makes the p- integral 
rapidly convergent (assuming that the p°-integral is carried out first). On the contrary, 
n B(|p°|) modifies the large-t behaviour of SV-^\t,r) significantly, since it is Bose-enhanced, 
n B(|p°|) ~ T/\p°\, for <C T. In fact, the contribution from this term does not vanish for 
t — > oo, but leads to a finite imaginary part for <5V> (oo,r) [ID] . 

In order to illustrate this behaviour, let us evaluate the term with (ly 1 ) numerically. 
An example is shown in Fig. [IJ We indeed observe that the imaginary part of the potential 
approaches a finite value at large t. 



4. Solution of the Schrodinger equation 



As argued in the previous section, the static potential in Eq. (|3.3[) should be evaluated in the 
limit t>r, yielding in dimensional regularization (cf. Eqs. (4.3), (4.4) of ref. [10] ) 



lim V^ 2 \t,r 



9 2 C F 

4-7T 



m D + 



exp(- 



ig 2 TC F 

4-7T 



(m D r) 



where the function 



dz z 



1 



sin(zx) 



zx 



(4.1) 



(4.2) 



/o {z 2 + l) 

is finite and strictly increasing, with the limiting values (p(0) = 0, (j)(oo) = 1. 

Before proceeding, it is appropriate to point out that by solving Eq. (|3.3p we only ac- 
count for a part of the 0(g 2 )-corrections, namely those which are temperature-dependent 
and change the i-dependence (or, after the Fourier-transform, the w-dependence) of the 
solution. Apart from these corrections, there are also other corrections, well-known from 
zero-temperature computations. In particular, the precise meaning of the mass parameter 
M should be specified; a matching computation between QCD and NRQCD p3] shows that 
it actually corresponds to a quark pole mass, whose relation to the commonly used MS 
mass is known up to 3-loop order [15] . Furthermore, the "normalization" of the NRQCD- 
representative of the electromagnetic current can be worked out by another matching com- 
putation: this relation is known up to 2-loop level [16]. In our language, this amounts to a 
radiative correction to the initial condition in Eq. (|2.9|) . Neither of these zero-temperature 
corrections plays an essential role for the thermal effects that we are interested in here, and 
consequently both will be ignored in the following. 
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4.1. General procedure 

Now, once Eq. (j3.3|) has been solved, we can extrapolate r — > 0, to obtain C>(f) = C>(i, 0). 
Symmetries indicate that C>(— i) = C>(i), whereby the Fourier transform from C>(t) to 
C>(lj) can be written as an integral over the positive half-axis. Recalling finally the relation 
of C > {uo) and the spectral function, Eq. (|2.2|) . we can write the latter as 



1 - e 



-/9w 



^°°dt {cos(wt)Re[c>(t,0)] - sin(wi) Im [(?> (i, 0)] } . (4.3) 



Concentrating on the non-relativistic regime, i.e. on frequencies close to the quarkonium 
mass, we write 

w = 2M + J , (4.4) 
with I a/ 1 <C M. It is also convenient to introduce 



<?>(*,r) = e- 



-»2Aft 



u(t,r) 



(4.5) 



Finally, we assume the point-split solution to be spherically symmetric (S-wave); in the 
following we denote it by u(t,r). Thereby Eq. (|4.3p becomes 



p(u>) 



1 - e 



-P(2M+uj') 



dt Uas(oj't)ReUp(t,0) — sin (a/ t) Im [-0(^0)1} 

^(t,0) = lim^l^ , 
and the Schrodinger equation reads 



where 



id t u(t,r) 



1 d 2 



+ vi 2) {oo,r) 



with the initial condition 



and the boundary condition 



Mdr 2 
u (0,r) = -6iV c r5 (3) (r) 
u{t,0) = . 



u(t,r) , 



(4.6) 
(4.7) 

(4.8) 

(4.9) 
(4.10) 



We note that the prefactor in Eq. (|4.6p can be set to unity, since we are in any case omitting 
exponentially small contributions ~ exp(— 2M/T). 



4.2. Discretised system 

In order to solve Eq. (|4.8p numerically, we discretise both the spatial coordinate r and the 
time coordinate We denote the spatial lattice spacing by a s and the temporal one by a t . 

1 Let us stress that this discretization is related to the solution of a classical partial differential equation; 
it has nothing to do with the regularization used in QCD. Indeed, Eq. (|4.ip assumes the use of dimensional 
regularization on the QCD side. 
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Furthermore, r max and t max are the maximal values of these coordinates; there are N s + 1 
spatial sites, and N t + 1 temporal sites, with r max = N s a s , i max = N t at- 

Let us start by discussing the discretization of the initial condition in Eq. (14, 9ft . In contin- 
uum, we can formally write 



oo _ f+1 

1 

d PP e ipr . (4.11) 



On the lattice, with r = na s , n = 0, 1, N s , a possible discretization of Eq. (|4.1ip . possessing 
formally the correct continuum limit at a s — * 0, is given by 



-*W(r) - -V f dp^sin(^) e ~ 

-(-l) n+1 . (4.12) 



t/c 

2 \ 2 n 



4n 2 - 1 



We will see in Sec. 14.31 from another angle that Eq. (|4.12[) indeed provides for a correct and 
very convenient discretization of the initial condition (once multiplied by — 6JV C ). 

As far as the spatial derivative in Eq. (|4.8p is concerned, we discretise it in the usual way: 

d 2 u(t,r) u(t, (n - l)a s ) - 2u(t, na s ) + u{t, (re + l)a s ) _ 

with the boundary condition in Eq. (|4,10p . Furthermore we also set the boundary condition 

u(t,N s a s ) = 0, (4.14) 

whose justification requires that we check the independence of the results on iV s (or r max ). 

As far as the discretization of the time derivative is concerned, the general issues arising 
are well described in §19.2 of ref. [T7j. Writing Eq. (|4.8p in the form 

id t u = Hu , (4.15) 

we use the "Crank-Nicolson method", which amounts to solving 

1 1 

(l + -iHa t )u{t + a u r) = (l - - iHa t )u(t, r) . (4.16) 

This method leads to an evolution which is accurate up to C(a 2 ), stable, and unitary (the 
last one provided that H were hermitean, which is not the case in our study). 

Given the solution for u(t,na s ), we then extrapolate for i()(t,Q) (cf. Eq. (|4.7p ) simply 
through 

lK*,0)ss^. (4.17) 
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4.3. Non-interacting limit in the discretised system 

The spectral function following from the discretization of Sec. 14. 2\ after the result has been 
inserted into Eq. (|4.6I) . can be found analytically in the free theory, if we take the limits 
o-t/o-s — ► 0, r max ,i max — ► oo. The solution is quite illuminating, so we briefly discuss it here. 
Let us start by introducing the notation 

2 fCL s p\ 1 
p= — sml — I, p=—sm(a s p). (4.18) 
a s v 2 > a s 

_____ (o) 

Then a general solution of Eq. (|4.8p [without V> (oo,r) and with the spatial derivative 
replaced by Eq. (|4.13|) ] can be written as 

u (t,r) = f W/aS ^ e-iP^/M+iprjr^ _ ( 41 g) 
J~n/a 3 2lT 

Satisfying the initial condition in Eqs. (|4.9p . (|4.12p requires 

^(p) = -6iV c -^ . (4.20) 

27T2 



V(t,0) = -6iV c -^ dppj3e-^ 2 '/ M , (4.21) 



Furthermore, extracting the function ijj(t,0) according to Eq. (|4.17p yields 

1 px/a 

-ir/a 

the Fourier-transform of which reads (cf. Eq. (|4,6p in the limit exp[— (2M + u')/T] = 0) 

QN C r . / , 4sin 2 (x/2) 



r / \ /a;\ ,/ , 4sin^(x/2)\ 
^ dxsin ( x)sin(-)^a^ -±±±) , (4.22) 



vra 2 . 



where uo' = lo — 2M. This integral can be carried out, with the outcome 

3N C M 2 _ . / 4 \ / \ ± 



= ^- 6>(u> - 2) {-^ + 2 - &) (u, - 2) . (4.23) 

We note that Eq. (|4.23p agrees exactly with Eq. (|2.7p . except that it is cut off sharply at 
(a) — 2) max = (2/a s M) 2 . For addressing the non-relativistic regime |u> — 2| <C 1 it is then 
sufficient to choose a s < 2/M for first estimates; at the end, one of course has to extrapolate 
a. s -> 0. 



5. Numerical results 

In a practical solution, we are not in the limit at/a s — * as in Sec. 14.31 but at is finite, 
and tmax) r max are finite as well. Then the time variable takes values t = nat, n = 0, ...,Nt, 
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while frequencies assume the values uj = 7rm/t max , m = —Nt,...,Nt- The Fourier-integral in 
Eq. (|4.6p is replaced by a discrete sum; to keep discretization errors at 0(a 2 ), we write it as 



/•tmax J r"i + "z 

J dtF(t) - - at Yl Finat) + Y,H 



Nt-l N t 



nat) 

n=0 n=l 



(5.1) 



For the parameter values needed we employ simple analytic expressions that can be ex- 
tracted from Ref. |18] . 

g2 " 91n(9.0 8 82T/A^) ' 31n(7.M7T/A^) ' ^N C = N { = 3. (5.2) 

We fix A^jg ~ 300 MeV; the width we will find does not depend significantly on this (see 
also Fig. 2 of ref. [10J). For the mass we insert the bottom quark mass, M ~ 4.25 GeV. We 
denote the "Bohr radius" by 

TB " ■ (5 ' 3) 

In the range of temperatures considered, g 2 Ci?/(47r) ~ 0.5. ..0.3, and r# ~ (4...6)/M. 

As typical values of the numerics-related parameters, we have used r max = 120 re, t max = 
r max , at = a s /5. The dependence on all of these parameters is beyond the visual resolution. 
By contrast, there is significant dependence on a s , given that discretization errors are only of 
order 0{a s ). We have consequently used several values and carried out a linear extrapolation 
to a s — > 0. A sufficient precision can be obtained, for instance, by using the values a s = re/12 
and a s = r#/24 for the extrapolation. 

The final result of our analysis is shown in Fig. [2l The curve "500 MeV <C T <C M" refers 
to the non- interacting result in Eq. (|2.7|) . 



6. Conclusions 

The purpose of this note has been to present a numerical estimate for the heavy quarkonium 
contribution to the spectral function of the electromagnetic current, based on Eqs. (|2.9j) . (13. 3ft . 
(|4.1|) , The conceptually new ingredient here is the inclusion of a thermal width through the 
imaginary part of the static potential in Eq. (|4.ip . 

The result we find, Fig. [21 shows a clear resonance peak which rapidly dissolves as the 
temperature is increased. Even though we do not expect the precise position and height 
of the peak to be quantitatively accurate, since higher-order perturbative corrections can 
be large in the temperature range considered (certainly up to 20%), it is comforting that a 
phenomenologically reasonable pattern arises from such a simple-minded computation. 

The result shown in Fig. [2] assumes that the spatial momentum of the electromagnetic 
current vanishes, q = 0. However, as discussed in Sec. 12. 1^ a non-zero q simply shifts the 
patterns horizontally by the center-of-mass energy q 2 /AM of the heavy quark-antiquark pair, 
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Figure 2: The bottomonium contribution to the spectral function of the electromagnetic current, 
divided by — M 2 , in the non-relativistic regime \u)/M — 2| <C 1. 



provided that q <C M. Furthermore, as also pointed out in Sec. 12. 1\ the dependence on the 
quark chemical potential is exponentially small in the range (M ± n)/T 3> 1. 

There has been a fair amount of interest in estimating the quarkonium spectral function 
from lattice QCD, mostly by making use of the so-called maximum entropy method [19j-|24j. 
Generically, these results show several resonance peaks, rather than just one as in Fig. [2j It 
has been suspected that the additional peaks may in fact be lattice artefacts. In spite of its 
own uncertainties, our computation seems to support such an interpretation. As far as the 
first peak is concerned, systematic uncertainties and different parametric choices do not allow 
for a quantitative comparison at the present time, but the patterns found on the lattice and 
in our study do appear to bear at least some qualitative resemblance to each other. 
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Appendix A. Auxiliary functions for Eq. (13.21) 

For completeness, we specify here the gluonic spectral functions that appear in Eq. (|3.2[) . In 
order to compactify the expressions somewhat, we introduce the notation 



y = 



P 



Then pt,Pe (cf. Appendix B of ref. [10] and references therein) can be written as 



Pt(p°,p) 

^t{v,p) 
(y 2 -l)p E (p°,p) 

&E{y,p) 

^E{y,p) 



9(y 2 - l)7rsign(y)5(A T (y,p)) + 



e(l-y 2 )T T (y,p) 



A 2 T (y,p)+T 2 T {y,p) 



2/„.2 



p (y 



m 



2 r 



2 , y 

V +2 



i-y 2 In 



7rm: 



%/(i-y 2 ), 



6'(y 2 -l)7rsign(y)<5(A B ( ? /,p)) + 



g(l-y 2 )r g (y,p) 
A|(y,p)+r|(y,p) 



2 , 2 

p + m D 



In 



y+l 



7rm; 



-y 



(A.l) 

(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 



It can be seen that there is in each case a contribution from the "plasmon" pole, as well as 
from the cut representing Landau damping. Restricting the integration to p° > thanks to 
reflection symmetry, the plasmon poles trivially yield 



dylC(y)5(A(y,p)) 



where yo > 1 is defined through A(y ,p) = 0, and 



/C(y ) 



\d y A(y ,p)\ 



\d y A T (y ,p)\ 
\d y A E (y ,p)\ 



yl 
y^- 
yo 

yo 



+ ^(i-2/o)in 



yl 



i 

l^yo + l 
2 yn-1 



yo - 1 



i ^ yo 

We note that the pole locations can be approximated as 



(A.8) 

(A.9) 
(A.10) 



yo 



mp 1 
P 



p » m D 
p <C m D 



for At , 



and 



1 + 2 exp 

mp 1 
P 



+ 1 



p 2> m D 
p <C m D 



for A s 



(A.ll) 



(A.12) 
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We finally remark that the integral over the angle between p and r in Eq. (|3.2p can be 
carried out, yielding 



+1 2 — e i P rz — e ~iprz 

dz = 2 

i 2 

-1 2 — g*?"" 2 _ e ~iprz 



sin(pr) 



pr 



(A.13) 



/+i 2 _ e v P rz _ g—iprz r , 

da: = = 2 cos(pr) — 1 + pr Sifpr) , (A.14) 

-i 2z z L J 

where Si(z) = Jq dt sm(t)/t. Then a two-dimensional integral is left over: the inner integra- 
tion over p°, the outer integration over p. 
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